Ising magnetism and ferroelectricity in Ca 3 CoMnOg 
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The origin of both the Ising chain magnetism and ferroelectricity in Ca3 CoMnOe is studied by ab 
initio electronic structure calculations and x-ray absorption spectroscopy. We find that CasCoMnOe 
has the alternate trigonal prismatic Co 2+ and octahedral Mn 4+ sites in the spin chain. Both the 
Co 2+ and Mn 4+ are in the high spin state. In addition, the Co 2+ has a huge orbital moment of 1.7 
fiB which is responsible for the significant Ising magnetism. The centrosymmetric crystal structure 
known so far is calculated to be unstable with respect to exchange striction in the experimentally 
observed Till antiferromagnetic structure for the Ising chain. The calculated inequivalence of the 
Co-Mn distances accounts for the ferroelectricity. 

PACS numbers: 71.20.-b, 78.70.Dm, 71.70.-d, 71.27.+a 



Among a variety of multiferroic materials discovered so 
far [3, Q, ferroelectric Ising chain magnet CaaCoMnOg is 
quite unique, because the ferroelectricity (FE) is driven 
by exchange striction in a collinear Ising spin chain con- 
sisting of the charge ordered transition-metal ions 
The spin chain has the alternate trigonal prismatic and 
octahedral sites 0, 0]. Special ft II antiferromagnetic 
(AF) structure is detected in CaaCoMnOe below Tn ~13 
K by neutron diffraction. However, the measured mag- 
netic moment of 0.66 ^b/Co and 1.93 /is/Mn is much 
smaller than the expected one of the normal high-spin 
(HS) Co 2 + (5=3/2) and Mn 4 + (5=3/2). This led Choi 
et al. to a conclusion that the Co 2+ is (surprisingly) in a 
low-spin (LS) state@. In contrast, the effective magnetic 
moment of fi e g =5.8-6.0 per formula unit (f.u.), ex- 
tracted from magnetic susceptibility measurements above 
T N @, suggests that both Co 2+ and Mn 4+ are in 
a HS state. Thus there is an apparent controversy be- 
tween those data, and the problem concerning the spin 
state of, in particular, Co 2+ ions seems to be still un- 
solved. Another important question is to understand the 
nature of the Ising magnetism and of the resulting ex- 
change striction, which are apparently crucial for the ap- 
pearance of FE in Ca 3 CoMnC>6. To this end, we carried 
out ab initio electronic structure calculations and x-ray 
absorption spectroscopy (XAS). We address the impor- 
tant issues including the Co/Mn site preference, their 
charge/spin/orbital states, the origin of the Ising mag- 
netism, and the exchange striction. 

Our ab initio calculations were performed by using 
the full-potential augmented plane waves plus local or- 
bital method (Wien2k code) [6(] . We took the experimen- 
tal centrosymmetric structure data of the rhombohedral 
lattice (i?-3c) which has in a hexagonal setting the lat- 
tice constant a=9.1314 A and c=10.5817 A [4, 7]. The 
calculations were done for different magnetic structures 



(TITI, Till, Till, and Tin orderings in the Co-Mn-Co- 
Mn chains). To study the exchange striction effect, we 
investigated the effect of internal atomic relaxation allow- 
ing the inversion symmetry to be broken in the till spin 
chain, as discussed below. The muffin-tin sphere radii are 
chosen to be 2.4, 2.1 and 1.4 Bohr for Ca, Co/Mn and O 
atoms, respectively. The cut-off energy of 16 Ryd is used 
for plane wave expansion, and 5x5x5 k-mesh for integra- 
tions over the Brillouin zone. To account for the strong 
electron correlations, GGA+U (the generalized gradient 
approximation Q plus Hubbard U) calculations Q were 
performed. J7=5.0 (4.0) eV for the Co (Mn) 3d electrons 
(with a common Hund exchange of 0.9 eV) are used .ac- 
cording to the calculations using a local-orbital basis [ifj . 
Note also that the choices of U value in the range of 2-7 
eV do not affect the conclusion made in this Letter. The 
spin-orbit coupling (SOC) turns out to be crucial and it 
is included by the second- variational method with scalar 
relativistic wave functions 

It is quite common that an octahedral Co 3+ ion is in 
a LS ground state. This may suggest that the small mo- 
ment of 0.66 ub /Co in CasCoMnOe could be due to the 
Co/Mn site disorder, i.e., an appreciable presence of the 
octahedral LS Co 3+ . To probe the Co/Mn site prefer- 
ence, we computed two structures either with the trigo- 
nal Co (Cotrig) and octahedral Mn (Mn oct ) or vice versa. 
The total energy results show that the former structure is 
more stable than the later by 0.33 eV/f.u. in GGA and 
more significantly, by 1.60 eV/f.u. in GGA+U+SOC. 
The energetically favorable structure has the HS trigo- 
nal Co 2+ and the HS octahedral Mn 4+ , while the un- 
favorable one indeed has the LS octahedral Co 3+ and 
the HS trigonal Mn 3 +. We plot in Fig. 1 the GGA 
and GGA+U+SOC calculated density of states (DOS) 
only for the favorable structure. The sharp DOS peak 
at the Fermi level, coming from the degenerate x 2 - 
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FIG. 1: (Color online) Density of states calculated by GGA 
(a) and GGA+U+SOC (b) for the tit! s P in structure of 
Ca3CoMn06 with the trigonal Co and octahedral Mn. It 
practically coincides with that for Till- The thin blue (bold 
red) curves refer to the majority (minority) spin. The Fermi 
level is set at zero energy. In (a) the Mn*J. has a closed t| 9 
shell centered at -1.5 eV, together with the e 9 bonding state 
around -5.5 eV. The Cotrig is in the high-spin 2+ state with 
the nearly degenerate 3z 2 -r 2 (do), x 2 -y 2 and xy (d±2) levels. 
In (b) the (x 2 -y 2 , xy) doublet splits due to the spin-orbit cou- 
pling, and the Hubbard U yields an insulating ground state 
with the minority-spin dnffe occupation. 

y 2 and xy (d±2) levels of the trigonal Co 2+ minority- 
spin d electrons, vanishes when going from the GGA to 
GGA+U+SOC solutions. This explains why, by open- 
ing a sizable gap of 2.4 eV, the GGA+U+SOC strongly 
favors the structure with the Co 2 r ~t" g and Mn^ t . The HS 
Co 2+ (Mn 4+ ) has a calculated spin moment of 2.64 (2.70) 
/is, both staying constant within 0.2 [is for U=2-7 eV. 

Having established the Co^g/Mn^ site preference 
from the above total-energy calculations, we turn to our 
XAS measurements to confirm it. The room temperature 
Co-L 2j 3 and Mn-L 2j 3 XAS of Ca 3 CoMn06 were collected 
at the ID08 beamline of the European Synchrotron Ra- 
diation Facility (ESRF) in Grenoble with a resolution of 
0.25 eV at Co-L 3 (at 780 eV). The sharp peak at 777.8 eV 
of the C0-L3 edge of single crystalline CoO and at 640 eV 
of the Mn-L3 of single crystalline MnO were used for en- 
ergy calibration. The spectra were recorded using the to- 
tal electron yield method by measuring the sample drain 
current in a chamber with a base pressure of 2 x 10~ 10 
mbar. Clean sample areas were obtained by cleaving the 
polycrystals in situ. 

Important is that XAS spectra are highly sensitive to 
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FIG. 2: The experimental and simulated x-ray absorption 
spectra of CasCoMnOe at the Co-Z/2,3 (upper panel) and the 
Mn-Z/2,3 (lower panel) edges, together with the experimental 
spectra of the references CoO, EUC0O3, MnO, and LaMnOs. 
CasCoMnOe turns out to have the high-spin trigonal Co 2+ 
and octahedral Mn 4+ . 

the valence state: an increase of the valence state of the 
transition metal ion by one causes a shift of the XAS 
£2,3 spectra by one or more eV towards higher ener- 
gies [ll], Ell- In Fig. 2 (upper panel) we see a shift 
of the center of gravity of the C0-L3 white line to higher 
photon energies by about 1.5 eV in going from the di- 
valent CoO to the trivalent EUC0O3. The energy po- 
sition of the Co line in CasCoMnOg is the same as in 
CoO, indicating a Co 2+ state. Also, a gradual shift of 
the center of gravity of the Mn-Z/3 white line to higher 
energies from MnO via LaMnOs to CasCoMnOe (lower 
panel in Fig. 2) evidences the increase of the Mn va- 
lence state from 2+ via 3+ to 4+. In fact the Mn-i^^ 
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edges of CasCoMnOg lie at the same energy position as 
in SrMn0 3 , LaMno.5Coo.5O3 [H, and LaMno.5Nio.5O3 
0, all of which have the Mn 4+ state. Thus, our XAS 
results confirm the Co 2+ -Mn 4+ state. 

To extract more detailed information concerning the 
charge and spin states from the Co-L2,3 and the Mn-L2,3 
XAS spectra, we have carried out simulations of the XAS 
spectra using the well-proven configuration-interaction 
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We studied a trigonal prism 
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cluster model 

C0O6 (an octahedral MnOg) cluster, including the full 
atomic multiplet theory and the local effects of the solid. 
Thus our model calculations account for the intra-atomic 
3d-3d and 2p-3d Coulomb interactions, the atomic 2p 
and 3c? spin-orbit couplings, the O 2p-Co 3d hybridiza- 
tion, and the proper local crystal-field parameters. The 
calculated Co-Z<2,3 XAS spectrum for the HS Co 2+ [curve 
(c) in Fig. 2] with the ionic trigonal crystal field inter- 
action A l 1 o tlic =0.75 eV reproduces the experimental spec- 
trum [curve (b)] very well [l7j]. In order to stabilize a LS 
state which was concluded by Choi et al [||, we would 
have to increase the A 1 ™ by nearly 3 times (2.8 eV), and 
the calculated spectrum for this case [curve (d) in Fig. 2] 
strongly disagrees with the experimental one [curve (b)]. 

In the lower panel of Fig. 2, one can see that the 
lineshapes of the experimental Mn spectrum [curve (f)] 
are well reproduced by the simulation [curve (g)] with a 
Mn 4+ (i| ff ) configuration in a local O^ symmetry. 

Thus, our ab initio calculations and XAS experi- 
ments have firmly established that the spin-chain mag- 
net Ca 3 CoMn0 6 



has the HS Co t 2 r j g and HS Mn 4 +, and 

3+ 



that there is no appreciable presence of the LS-Co J./HS- 
Mn^. Obviously, the magnetic moment of 0.66 /is/Co 
measured by neutron diffraction Q below Tn ~13 K is 
much smaller than our theoretical value [the total cal- 
culated moment being about 4.3 ^ts/Co with the spin 
(orbital) contribution of 2.6 (1.7) hb, see more below] 
and than the value of fj, e s=5. 8-6.0 /is obtained from 
the high-temperature magnetic susceptibility 0, [B[ . This 
may be partially due to strong fluctuations in quasi one- 
dimensional chains, with frustrated interchain interac- 
tions. Another important factor may be a close prox- 
imity of different types of magnetic orderings, see below. 
This question deserves further study. 

As seen in Fig. 1(a), the trigonal Co 2+ has the nearly 
degenerate 3^ 2 -r 2 and (x 2 -y 2 , xy) levels. Due to the 
in-plane character of both the x -y and xy orbitals, 
their strong Coulomb repulsion prevents their double oc- 
cupation in the minority-spin channel. Therefore, the 
minority-spin 3z 2 -r 2 orbital is fully occupied and the 
minority-spin x 2 -y 2 and xy are half-filled for the HS 
Co 2+ ions. Due to the quasi one-dimensional character 
along the c-axis chain, a naive x 2 -y 2 /xy planar orbital 
ordering does not gain any energy (compared with either 
an x 2 -y 2 1 x 2 -y 2 or xy/xy orbital ordering), as proved by 
our ab initio calculations. In contrast, an efficient way to 
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FIG. 3: Relative energy as function of the Co displacements 
calculated by GGA+U+SOC for the experimental lattice. 
Triangles stand for the data points, which are fitted into a 
parabolic curve. The energy minimum shows the Co displace- 
ment of 0.006 A, which leads to alternate short-long-short- 
long Co-Mn distances (the difference being 0.012 A) in the 
Cot-Mn|-Co|-Mn| chain as depicted schematically in inset. 

gain the energy is SOC. When the SOC is included, the 
(x 2 -y 2 , xy) doublet splits into d 2 and d_2, and the gain 
of the full SOC energy is calculated to be about 70 meV. 
As a result, a huge orbital magnetic moment of 1.7 /is is 
generated at the HS Co 2 ^~ g sites with the minority-spin 
doO?2 occupation [Fig. 1(b)], and the SOC firmly fixes the 
total magnetization (with parallel spin and orbital con- 
tributions) along the c-axis chain direction. Note that 
the orbital moment stays constant within 0.1 ub in our 
GGA+U+SOC calculations for [7=2-7 eV. Therefore, the 
peculiar trigonal crystal field and the SOC are responsi- 
ble for the significant Ising character of CasCoMnOfi; cf. 
similar situation in the isostructural CasCoRhOg [191 ]. 

Having established the picture about the Co/Mn site 
preference, their charge/spin/orbital states and the sig- 
nificant Ising magnetism, we turn now to the discussion 
on the intra-chain magnetism and the exchange striction 
leading to ferroelectricity. Since a centrosymmetric struc- 
ture of CasCoMnOe does not carry ferroelectricity, we 
now look whether the symmetry can be lowered by allow- 
ing for atomic displacements within the experimentally 
observed [1] |t|| magnetic structure of the Co-Mn-Co- 
Mn Ising spin chain: the Co-Mn bonds would be inequiv- 
alent in such a non-centrosymmetric structure leading to 
chain dimerization, and this will give rise to ferroelec- 
tricity. We calculated using GGA+U+SOC the total en- 
ergy of the system by firstly shifting the Co ions along 
the chain making Co-Mn distances unequal (alternating) 
[20I ]. The results are shown in Fig. 3. We see that indeed 
the lattice in TTII structure would relax to a state with 
alternating Co-Mn distances (the difference being 0.012 
A), which means the appearance of FE in the system, cf. 
HHH HE Hi- Th e calculated spin and orbital moments 
remain unchanged. Furthermore, a triclmic lattice (PI) 
was tested to study a likely breaking of the 3-fold rotation 
symmetry. A structural relaxation using GGA+U+SOC 
(with the Co ions fixed at the above optimized position) 
confirms the small FE distortion of the Co-Mn chains and 
gives also a small distortion of the Co-0 (Mn-O) bonds 
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[only about 0.01 A out of 2.15 A (1.92 A) in bondlength] . 
We would like to point out that those small distortions 
are not inconsistent with the experiment Q since they 
are so small that they would require a special effort for 
their detection. 

We note that in the GGA+U+SOC calculations the 
ft II an d Tit I magnetic structures for the Co-Mn-Co- 
Mn spin chain are the two lowest-energy states, with the 
former being higher in energy. The energy difference is 15 
meV/f.u. for the experimentally detected crystal struc- 
ture Q and is reduced to 7 meV/f.u. for the relaxed 
triclinic lattice (PI). At the moment it is not clear how 
to stabilize the experimentally observed ft II magnetic 
structure 0] in the GGA+U+SOC calculations. Since 
the energy difference is rather small and the intrachain 
exchange interactions rather weak, we speculate that the 
couplings between the chains may play an important role. 
The present calculations assumed a simple ferromagnetic 
interchain interaction, but it is known that there are 
AF interchain couplings in the isostructural Ca3Co206 
[25| and CasCoRhOe compounds with strengths not 
much smaller than the intrachain ones. These finite AF 
interchain couplings bring about magnetic frustration, 
and one needs to investigate a large set of intra- and 
interchain magnetic structures in order to explain theo- 
retically the magnetic ground state of CasCoMnOg- 

In summary, we find the site preference of the high- 
spin trigonal Co 2+ and the octahedral Mn 4+ in the Ising 
chain magnet CaaCoMnOg, using the ab initio calcula- 
tions and x-ray absorption spectroscopy. The Co 2+ has 
the very stable minority-spin <io^2 occupation due to the 
peculiar trigonal crystal field, inter-orbital Coulomb re- 
pulsion and the spin-orbit coupling. Thus, a huge orbital 
moment of 1.7 hb is predicted and the significant Ising 
magnetism is well accounted for. Moreover, our calcu- 
lations indeed found that a structural relaxation due to 
the exchange striction decreases the energy of the exper- 
imentally observed till magnetic ordering and leads to 
the observed ferroelectricity. However, the mechanism of 
the full stabilization of the Till structure as the ground 
state calls for a further study. 
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